Phase Coherence in a Driven Double- Well System 
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We analyze the dynamics of the molecular field incoherently pumped by the photoassociation 
of fermionic atoms and coupled by quantum tunnelling in a double-well potential. The relative 
phase distribution of the molecular modes in each well and their phase coherence are shown to build 
up owing to quantum mechanical fluctuations starting from the vacuum state. We identify three 
qualitatively different steady-state phase distributions, depending on the ratio of the molecule- 
molecule interaction strength to interwell tunnelling, and examine the crossover from a phase- 
coherent regime to a phase-incoherent regime as this ratio increases. 
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The theoretical analysis of interacting quantum fields 
normally requires the introduction of some approxima- 
tion schemes that limits the detail with which they can 
be understood. In the context of many-body theory, 
these approximations often take the form of a factor- 
ization (and possibly linearization) of the higher order 
correlations for the field operators such as in mean-field 
theory y| . There are of course exceptions to this state of 
affairs, many of them in the realm of quantum optics Q. 
One such example is the micromaser, which is amenable 
to an almost exact quantum description that permits a 
detailed dynamical understanding of the Maxwell field 
inside a high-Q microwave resonator 0, • 

Ultracold matter-wave fields with adjustable interac- 
tions [5j are now under remarkable experimental control 
and are becoming more widely available. This makes 
them an ideal test bed for detailed dynamical studies 
of interacting quantum field theories. In this letter, we 
discuss one such example, a matter-wave analog of cou- 
pled micromasers that can be realized using experimental 
techniques recently developed in the quest for molecular 
Bose-Einstein condensation The dynamics of that 
system is amenable to an essentially exact analysis that 
shows in particular the build-up of phase coherence be- 
tween the coupled molecular states. As such, it sheds new 
light on the dynamics of quantum fields in a system that 
should soon be amenable to experimental realization. 

A major current thrust in atomic, molecular and opti- 
cal science is the study of molecular Bose-Einstein con- 
densates (BEC) formed via either photoassociation 0, 
or a Feshbach resonance 0, 0, 0|. These efforts have 
recently culminated in the realization of molecular con- 
densates of 40 K 2 and 6 Li 2 molecules 0- 

The formation of molecules by photoassociation in op- 
tical lattices has been studied theoretically by several au- 
thors ■ The large energy separations that are possible 
between the lowest and second Bloch band of the lattice 
allow one to restrict the center-of-mass states of atoms 
and molecules to the lowest energy Wannier state of each 
site, thereby avoiding many of the difficulties associated 
with free space Search et al. 0] recently showed 



that the photoassociation of fermionic atoms into bosonic 
molecules in a lattice in the Mott-insulator regime (l4Hl5| 
can be mapped on to the system of the micromaser, a de- 
vice that exhibits a number of non-classical features Q • 

The present letter extends this model to include inter- 
well tunnelling in a double-well system 0, 0, 0] and 
analyzes the full dynamics of the molecular field. We 
examine the build-up of the relative phase between the 
localized molecular states of the wells due to the com- 
bined effect of inter-well tunnelling and the incoherent 
addition of molecules from photoassociation. 

We identify three regimes, characterized by different 
orders of magnitude of the ratio of the two-body colli- 
sion strength to the inter- well tunnelling coupling |lfij . 
For small values of this ratio the relative phase of the 
two molecular modes becomes phase locked at a fixed 
value, while for larger values that phase remains random. 
This crossover of the non-equilibrium steady state from 
a phase coherent regime to the phase incoherent regime 
is reminiscent of the superfluid-Mott insulator transition 
for the ground state in an infinite lattice flil llfij . How- 
ever, in contrast to that system, this behavior occurs now 
in an open system, with incoherent pumping and damp- 
ing of the matter-wave field. 

We consider the two-photon stimulated Raman pho- 
toassociation of fermionic atoms of mass rrif and spin 
(er = 1, 2) trapped in a double- well potential into bosonic 
molecules of mass m& = 2m / 0, 13 ■ We assume that the 
system is at zero temperature and that the number of 
fcrmions of each component that occupy each well is less 
than or equal to 1 at all times. The fcrmions, therefore, 
occupy only the lowest energy level. The molecules are 
confined in the same potential and occupy only the lowest 
energy level of each well. 

The effective Hamiltonian for the atom-molecule sys- 
tem in the lowest energy level of the wells is H = 
Z)i=j,r(-^oi + Hn) + H T , where H 0i = h(uj b + 5)nj + 
fiuf(h li + h 2i ) + ^hUbhitfit - 1) + hU x hi(hii + n 2i ) + 
h,Ufhuh 2i , Hn = Tix{t)b\c u c 2i + h.c, and H T = 



-hJbbjbr 



hJf(c u ci r + c 2l c 2r ) + h.c. Here c„i and bi are 
the annihilation operators of fermionic atoms and bosonic 
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molecules in the left or right wells, i = l,r, respec- 
tively. The corresponding number operators hi = S^S; 
and fieri = c'oiCoi have eigenvalues n, and n a %, respec- 
tively, and u>b and ojf are the energies of the molecules 
and atoms in the isolated wells (Jb,f = 0). The terms 
proportional to Ub, U x , and Uf in Hqi describe on- 
site two-body interactions between molecules, atoms and 
molecules, and atoms, respectively. The parameters Jb 
and Jf in the tunnelling Hamiltonian, Ht, are the sin- 
gle molecule and atom tunnelling rates. The interac- 
tion Hamiltonian, Hu, describes the conversion of atoms 
into ground-state molecules via two-photon stimulated 
Raman photoassociation. The photoassociation coupling 
constant is proportional to the far off-resonant two- 
photon Rabi frequency associated with two nearly co- 
propagating lasers with frequencies u)\ and u>2 |8j, and 8 
is the two-photon detuning between the lasers and the 
energy difference of the atom pairs and molecules. 

The dynamics of the molecular field in the double-well 
system is governed by four mechanisms: (1) The injection 
of pairs of fermionic atoms inside the double- well during 
time interval T when X = 0, after which the atoms in each 
well are in the state |e,-) = c^Cj^O) with unit probability 
[l3T | ; (2) The unitary time evolution of the molecular field 
subject to the Hamiltonian Hb during the time interval 
T, with 



Hb = —hJb(bjb r 



blk) + h(U b /4)(hi - h r 



(1) 



In Eq. QJ, we have neglected terms that are functions 
only of hi + h r . This is justified as long as the ini- 
tial density matrix is diagonal in the total number of 
molecules in the two wells (see below); (3) The molecu- 
lar damping at rate 7 during the intervals T |13I |. mod- 
elled by a master equation with a Liouvillian Cipb = 
-(^/2)[b\biPb ~ 2bipbb\ + pbb\bi] for each well, p b being 
the reduced density operator for the molecules; (4) The 
switching on of the photoassociation lasers in a train of 
short pulses of duration r and period T + r. We as- 
sume that r is much shorter than all other time scales 
in this model, r <g; Jj^^j^ 1 , so that damping and tun- 
nelling may be ignored while the photoassociation fields 
are on. The atom-molecule conversion is described by 
Fi(T)p b = Tr a [Ui(T)pab(t)UI(r)],wh.eie p ab is the total 
density operator and Tr a [ ] denotes the trace over the 
atomic variables, Ui(r) — exp (—ihir/h) being the evo- 
lution operator for the Jaynes-Cummings-like Hamilto- 
nian of each well. It describes the interaction between 
the molecular field and a fictitious two-state system for 
fermionic atoms |e») and \gi) = |0) fl^j . 

hi = h {uj b + U x ) hi + h (ujf + U x hi) a zi 

+ h (x{t)b\a-i + x*(t)6iff+<) + ^U b hi(hi ~ 1) (2) 

where = = \gi)(ei\ and a zl = |e,)(ei| - \gi){gi\. 
In Eq. (|2J) we have dropped constant terms and made 



the redefinitions ujb + 8 — > Ub and Uf + Uf/2 — > ujf. 
For x — const, Eq. J5J) can be solved within the 
two-state manifolds of each well {|ej, rij), |<7,-, + 1)}. 
Within each manifold, the resulting dynamics is in the 
form of quantized Rabi oscillations at frequency 7i ni = 



[2u f ~Lu b + (2U X - U b ) n,X + 4| X | 2 K + 1). 
The photoassociation of atoms into molecules during 
the intervals r leads to the build-up of the molecular 
field. The coarse-grained master equation for the reduced 
molecular density operator, valid for T 3> r, is given by 



Pb 



E fa 

j=l,r 



■T-HW-Ij}] Pb-^ 



Hb, Pb 



(3) 



The density operator for the molecular field in the num- 
ber representation for the left and right wells is given 
by Pb = En,,n r ,m i ,m r / 5 ( ri i' n '-;m i! m r )|ni)|n r )(m r |(rni|. 
The initial condition for the molecules is taken to be 
the vacuum state. Because the molecular pumping and 
decay is the same in both wells, p remains diagonal in 
the total number of molecules in the two wells, N = 
m + n r = mi + m r , for all times. Consequently, the 
only nonzero terms in p(ni, n r ; mi, m r ) are those with 
ni + n r = mi + m r . 

The master equation J3J depends on six independent 
parameters: the pump parameter, = \J N cx \x\t; the 
number of photoassociation cycles per lifetime of the 
molecule, N ex = I/7T; the two-body collision strength 
and tunneling coupling strength per decay rate, u b = 
Ub/') and tj = Jb/l', and finally, the (non-)linear de- 
tuning parameters, 77 = (2a;/ — !Jb)/2|x| (j3 = (2U X — 
Ub)/2\x\) from lZ ni . In the following, we consider for 
simplicity only exact resonance, r\ = [3 = 0, and fixed 
N ex = 10. The results presented below do not depend 
in any qualitative manner on the specific value of these 
parameters. 

Introducing the angular momentum representation 



J. T, 



= bjb 



J+ = JL = 

and J 2 = (N/2)(N/2 + 1) 
is the total number operator 

({->* = E nil „ r vW 

(Jz) = En^nJ"! " n r )p{ni,n r ;ni,n r )/2 



J z = (blk - b\b r )/2, 
where N = hi + h r 
we have that (</+) = 
l)n r p(ni,n r ;ni + l,n r — 1), 

Since the ini- 



tial state of the molecules in each well is the same and 
Hb is invariant with respect to the interchange I <-> r, 
the density matrix is invariant with respect to ni «-> n r 
and mi <-> m r . As a result, the reduced density matrices 
for the left and right wells are identical. This leads to 
the same single- well molecule statistics for the two wells. 
The symmetry of the density matrix with respect to the 



two wells furthermore implies that (J z ) = (J y 



0. 



Figure HJa) shows the normalized steady-state first- 
order coherence (J x ) / (hj) as a function of Ub/tj for 
= 7r and tj = 2.5. (J x ) is suppressed in both the weak 
and the strong two-body coupling limits, and has an ex- 
ternum at \ub\/tj r*> 0.6. In the strong coupling regime, 
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FIG. 1: (a) (Jx)/{nj) and (b) gjy versus ut/tj for 6 = tv and 
tj - 2.5. 

\ub\/tj ^S> (N) ~ 10, the nonlinearity in H b dominates 
and reduces the coherence between the localized states 
of each well. We note that the average occupation num- 
bers for each well are relatively unaffected by Ub/tj, with 
(hj) = (N) /2 = 4.78 - 4.87 for \u b \/tj = 10 2 - 1CT 2 - 5 . 

The origin of the mutual coherence between the two 
molecular modes can be determined from the equation of 
motion for (J x ), 

(Jx) = ^ y/(ni + l)n r p R (ni,n r ;ni + l,n r -l) 

= -U b (jyJ Z + J Z Jy) - 7(J X ) 

+ 2T _1 ^2 Vn~ip ni p R (ni,n r ;ni + l,n r - 1), (4) 

m ,n r 

where p ni = - 1) + VhJ+2S% +1 

with C£ = cos (R n r / 2) cos (R m r/ 2), = 
sin (R n r/2) sin (R m r/2), and pn is the real part of 
p. In deriving Eq. (J3J we used the I r symmetry of p 
and p(ni,n r ;mi,m r ) = p*(mi,m r ;ni,n r ). 

In the absence of collisions, C4 = 0, it can be shown 
from Eq. © that pii(ni,n r ;ni + q,n r — q) with q odd 
couple only to themselves and to the imaginary part 
of p(ni,n r ;ni + p,n r — p) with p even. For wells in- 
coherently pumped at equal rates, this implies that 
Pft(rii,n r ;ni + l,n r — 1) = for all times. Alterna- 
tively this can be understood by noting that the expec- 
tation value (J x ) corresponds to the difference in occu- 
pation numbers between the in-phase, b s = (b\ + 6 r )/v2, 
and out-of-phase, b a — (b\ — b r )/y/2, states of the local- 
ized states of each well, J x = b\b s — b' a b a . These states 
are equally populated since the bandwidth of the pho- 
toassociation pulse is larger than their energy splitting, 
1/t> Jfc. Hence (J x ) — for U b — 0, and the creation 
of cross-coherence between the two wells is due solely 
to two-body collisions. We also remark that a semiclas- 
sical treatment in which (J y J z + JzJy) is factorized as 
2(J y )(J z ) [18( results in (J x ) = for all times and all 
values of u b /tj. Hence, the build-up of ( J x ) is a purely 
quantum-mechanical effect due to quantum fluctuations. 

The pump and decay mechanisms act identically on 



the in-phase and out-of-phase states, hence the fact that 
(J x ) has the sign of U b results from the unitary time 
evolution from H h , which gives (J x (t)} = (u b /2tj)(J^(t)), 
if (J x (0)) = (J 2 (0)) = 13. Since (J 2 ) > 0, it follows 
that the sign of (J x ) is determined by that of U b . 

Figure njb) shows the steady-state second-order cor- 
relation function, g£' = (hih r ) / (hi) (h r ) as a function 
of \u b \/tj for 6 = 7r and tj = 2.5. It is related to 
the variance of the relative number difference, such that 
(AJ Z ) 2 = (Am/2) 2 - (hi)(h r )(g { ^ - l)/2 + (An r /2) 2 . 
In the strong two-body coupling limit, the molecular 
states are uncorrelated gj® — 1, and (AJ Z ) 2 becomes 
a sum of variances for independent localized states of 
the molecules. In the weak coupling region, \u b \/tj <C 
(N)^ 1 , the two molecular states are anti-correlated, 
gj^ < 1, and the variance (AJ Z ) 2 is greater than that of 
uncorrelated states. The enhancement of relative num- 
ber fluctuations indicates the locking of the relative phase 
between the two wells. 

To investigate the relative phase distribution of the 
molecular field in the two wells, we consider the differ- 
ence of the two single-mode Pegg-Barnett phase opera- 
tors Specifically, we introduce the (s + 1) 2 orthonor- 
mal phase states 

s s 

mdr) = JVTT) £ £ e-'V^lnOK), (5) 

where 0; = 2nl/(s + 1) and 9 r = 2irr/(s + 1) (l,r = 
0,l,---,s) are phase variables and s is a finite value. 
Projection onto these phase eigenstates gives the quan- 
tum phase distribution for the two wells, P(8i,8 r ) = 
Tr{\6i,6 r )(6i,6 r \pb}. Since the density matrix is diago- 
nal in the total number of molecules, rt; + n r = mi + m r , 
the phase distribution depends only on the relative phase, 
4>l,r = 8i — 6 r , which has a width of 47r in the phase co- 
ordinate. In order to obtain a mod(27r) distribution for 
the relative phase, we average over the total phase coor- 
dinate, 9i + 9 r , and map the mod(47r) distribution into a 
mod(27r) one following the method of Ref. |21| . Using a 
redefined relative phase variable, <f> n = 2itn/(s + 1) — tt 
(n = 0, 1, • • • , s) gives the mod(27r) distribution 

p ( ( l ) n) = —-r V V e _lfc ' 6 "p(n;,n r .;n / +fc,?i r -fc). 
s + 1 * — ' L — ' 

ri( ,n r — k=—ni 

Figure |21 shows the evolution of the relative phase dis- 
tribution P(4> n ) in three different regimes (a) u b /tj — 
0.0032, (b) Ub/t,, = 0.5623, (c) u b /tj = 56.23, for 6 = tt, 
tj = 2.5, and s = 40. Since the vacuum state is taken 
as the initial state, the relative phase at t = is ran- 
domly distributed, P{<t> n ) = l/(s + l). For \u b \/tj < (N), 
corresponding to Figs. Efa) and (b), P((f> n ) develops a 
peak around and/or ±7r in the characteristic time 
needed to reach a steady state Q. As previously dis- 
cussed, for weak two-body interaction the in-phase and 
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FIG. 2: Time evolution of P(<j> n ) for 9 = n, tj = 2.5 and for 
(a) u b /tj = 0.0032, (b) u b /tj = 0.5623, (c) u b /tj = 56.23; 
(d) (Acf>) 2 of the steady state. 



out-of-phase modes have equal population. This leads to 
the bimodal phase distribution with peaks around both 
and ±7r. For moderate Ub/tj, the relative phase locks 
in time around (±7r), for repulsive (attractive) two- 
body interactions, see Fig. |2fb). In contrast to these 
two regimes, when \ub\/tj ^> (N), the relative phase 
distribution becomes almost random, and the localized 
modes in the two wells evolve independently of each 
other. The steady-state phase variance, (A0) 2 , is shown 
in Fig. Hfd) as a function of the ratio Ub/tj. Consistently 
with Fig. Ufa-c), it exhibits large fluctuations due to the 
bimodal phase distribution in the case of Ub/tj < (N)^ 1 , 
and becomes narrow in the region (N)^ 1 < Ub/tj < (N). 
When Ub/tj > (N), the variance approaches the value 
(A0) 2 = 7r 2 /3 corresponding to a uniformly distributed 
phase (s — * oo). 

The three regimes of phase distributions correspond 
to different orders of magnitude of the ratio Ub/tj, which 
characterizes the behavior of the double- well system . 
The crossover of the non-equilibrium steady state from a 
phase-coherent regime to the random-phase situation is 
reminiscent of the superfluid-Mott insulator phase tran- 
sition for the ground state of an optical lattice [Til fl5| . 
Since we consider just two sites, however, there is no 
sharp transition between these regimes. 

In this paper, we have studied the quantum dynamics 
of the bosonic molecules in a double-well system in the 
presence of both incoherent pumping and damping. We 
have shown that the phase coherence of initially indepen- 
dent molecular states builds up due to quantum mechan- 
ical fluctuations. We identified three qualitatively dif- 



ferent phase distributions for the non-equilibrium steady 
states. The crossover from the phase coherent regime 
to phase incoherent regime occurs as the ratio \ub\/tj 
increases, similar to the superfluid-Mott insulator phase 
transition in an optical lattice. 
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